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Abstract 

The low temperature and large volume effects in the d = 2 + 1 antifer- 
romagnetic quantum Heisenberg model are dominated by magnon exci- 
tations. The leading and next-to-leading corrections are fully controlled 
by three physical constants, the spin stiffness, the spin wave velocity and 
the staggered magnetization. Among others, the free energy, the ground 
state energy, the low lying excitations, staggered magnetization, stag- 
gered and uniform susceptibilities are studied here. The special limits of 
very low temperature and infinite volume are considered also. 
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1 Introduction 



The low energy properties and the leading finite size and temperature effects of the 
d = 2 + 1 quantum anti-ferromagnetic (AF) Heisenberg model 



H = JJ2SnS^^^, si = 3(3+1) (1.1) 

n./j, 

are governed by magnon excitations. The interaction of massless excitations whose 
existence is due to spontaneous symmetry breaking [Q], is strongly constrained by 
symmetry principles 0. In particle physics it was the interaction of soft pions where 
many of the general features of this dynamics were first observed |]^, ^, These 
investigations received new impetus recently by the application of the powerful 
method of effective Lagrangeans 0, ^] which provides a systematical way to calcu- 
late higher order corrections (chiral perturbation theory). Gasser and Leutwyler and 
their followers treated not only low energy Goldstone boson physics but also finite 
temperature and finite size effects in QCD 10|. This technique can be imme- 



diately generalized to other quantum field theories and critical statistical systems 



TT| and provides a systematic alternative to other related methods [|T^, [T^, |T^ . 



The special, constrained feature of magnon interactions in the AF Heisenberg 
model has been emphasized and different leading order results were obtained in sev- 



eral recent publications [|15| |T6|, ^7j. Using current algebra techniques or intuitive 
reasoning it would be very difficult to go beyond the leading order. On the other 
hand, this problem becomes a relatively simple book-keeping if the technique of 
chiral perturbation theory is used. The power of chiral perturbation theory was 
illustrated in this context in ref. |18[ by calculating the first correction to the corre- 
lation length at low temperatures. In this paper we present further results beyond 
leading order. The temperature and size dependence of the free energy density (in- 
ternal energy density, specific heat), the size dependence of the ground state energy 
density, the energy of the low lying excitations and its size dependence will be inves- 
tigated. We calculate the staggered magnetic field dependent part of the free energy 
density (magnetization and susceptibility) and study the behaviour of the uniform 
susceptibility. Additionally, we discuss the zero temperature and the infinite volume 
limits which reveal special physical properties and pose new technical problems. 

The results of chiral perturbation theory refiect the symmetries of the underlying 
model only and depend on free physical parameters whose number is growing with 
the order of the calculation. As it is well known, the leading order results depend 
on three physical parameters: the spin stiffness constant ps, the spin wave velocity 
c and the staggered magnetization Ais- It is less well known, perhaps, that for the 
quantities we study here the next-to-leading corrections can be expressed in terms 
of the same constants, no new parameters enter (see ref. |11| and Section 4). As we 



shall discuss, these predictions are free of cut-off effects also. 
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The three parameters p^, c and Aig are fixed by the underlying model as the 
function of the coupling J. An effective way to fix these constants is to compare 
the results of chiral perturbation theory with the corresponding numerical results 
using the Hamiltonian operator in eq. This method has been applied earlier 

for classical ferromagnetic models and field theories ||T9| . 

The classical 0(3) ferromagnetic Heisenberg model in li = 3 in the broken phase 
is described also by the results of this paper. In this case we have to put c = 1 and 
interpret Lt = and Lg as giving the size of the euclidean box Lt x Lg x Lg on 
which the model is defined. Actually, we give all the results for the symmetry group 
O(A^), = 3 corresponding to the model in eq. ( |1 . 1| ) . 

The paper is organized as follows. For the readers convenience we summarize 
the results in Section 2. Concerning the notations, results and region of validity, 
this section is self-contained. Some of the results discussed in Section 2 occurred 
earlier in the literature in the context of classical ferromagnets or field theories. In 
this case our role is only to insert the proper spin wave velocity dependence and to 
raise attention to the existence of these results. In Section 3 we try to clarify some 
problematic results which occur in the literature. In Section 4 we construct the 
effective Lagrangean which is needed to derive the free energy density / = /(T, Lg) 
up to next-to-leading order. The effect of a finite cut-off is discussed in some detail 
there also. In Section 5 we calculate f(T, Lg) up to next-to-leading order. Although 
this calculation is a priori valid for hc/{TLg) =0(1) only, the — >■ oo limit gives 
correctly the temperature dependence of the infinite volume free energy density at 
low temperatures. The formal T ^ limit, however, does not lead to the correct 
volume dependence of the ground state energy, as we discuss in Section 5. The vol- 
ume dependence of the ground state energy, the lowest lying excitations and their 
finite volume corrections are considered in Section 6. In this section we regularize 
the effective Lagrangean on a lattice which opens the possibility to discuss the issue 
of power divergences in chiral perturbation theory. We show that — as expected — 
all these divergences can be absorbed through the standard renormalization proce- 
dure and no mixing between infrared and ultraviolet problems occurs. In Section 7 
the dependence of the free energy density on a weak staggered magnetic field, the 
magnetization and susceptibility are considered. Since this problem has been treated 
before we can keep this part brief. Very low temperatures require special care again 
as discussed in Section 8, where the volume dependence of the staggered susceptibil- 
ity at T = is calculated. In Section 9 the uniform susceptibility is discussed. In the 
Appendix A we consider specific lattice momentum sums and integrals. Appendix 
B contains some asymptotic expressions for the shape coefficients. 

Chiral perturbation theory is an exact low energy expansion. It summarizes the 
symmetry properties of the model in eq. ( |1 . ![ ) in a concised way and it is valid for any 
spin S. Consequently, if we replace in any quantity considered here the parameters 
Ps, c and A4g by their spin wave (1/5*) expansion, then the spin wave expansion of 
this quantity is obtained. The relation between these two expansions is discussed 
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further in Section 3. 

We made an effort to stay close to the notations used earher in the hterature. 
Since the model in eq. ( |1.1| ) gave the motivation for this work, we denoted by ps 
(rather than by F^, or T) the leading non-magnon scale. We used the index s both 
to denote staggered quantities and to refer to 'spatial', like spatial volume and size. 
Certain parts of the calculation are carried out for general dimension d, in which 
case we kept d as a parameter in the formulas. In the calculations we used the 
convention = c = 1" and restored the spin wave velocity dependence only at the 
end. At some places we introduced more convenient notations with respect to the 
preprint version of this paper (BUTP-92/46). 



2 Results 



Let us introduce a staggered external field = (-B^, 0, 0), and a uniform magnetic 
field0 hu = 0, 0) and generalize eq. (pTTP to 

H = JJ2SA+^ - B, ^(-l)ll"lls„ - K E Sn . (2.1) 

n,fi n n 

We shall define the field H, as 



H. = ^ , (2.2) 



where a is the lattice unit of the square lattice over which eq. (|2.1|) is defined. 
We shall consider a system of spatial size Ls x Ls at temperature T. The spin- 
stiffness, the spin wave velocity and the staggered magnetization of the infinitely 
large system at zero temperature and at zero external fields are denoted by ps, 
c and A4s, respectively. These quantities are fixed by the Hamilton operator in 



eq. (|2.1| ). On dimensional grounds ps and he/ a are proportional to J. It will be 



useful to introduce the notation 

In the path integral formulation we shall work in a euclidean box Lt x Lg x Ls, 
where Lf = hc/T. We shall refer to the / = 0(1) and / ^ 1 cases as "cubic" and 
"cylinder" geometry, respectively. Concerning dimensions, we have 

erg 

dim J = erg, dim Hs = — 5-, dim Lg = cm, dim T = erg, 

cm"' 

dim hu = erg, dim ps = erg, dim M.s = 1, dim 1 = 1. (2.4) 



The magnetic moment is included in the definition of 
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The free energy density / is defined as 

Z = exp (^-^/^ , (2.5) 

where Z is the partition function. As explained in the introduction, the following 
results apply to classical 0(iV) symmetric ferromagnetic systems also, while eq. 
corresponds to = 3. In general, magnons dominate the physical properties of the 
system at low energies, large volumes, low temperatures and weak external fields. 
More detailed conditions will be given together with the results. 



2.1 i^g— dependence of the basic parameters (Lg = oo, T = 0, 

^0,K = 0) 



In the presence of a nonzero staggered field Hg the magnons (they have a relativistic 
dispersion relation) pick up a mass . Let us introduce the notation M: 

i^^Mf=^M^Hg. (2.6) 



Ps 



We shall consider weak external staggered fields such that 



(2.7) 



Denoting the staggered magnetization and the spin stiffness in the presence of the 
field Hs^^hy and , we have 



M 



H 



H 



P. 



Ps 

M 



^ N-l c^M ^.^^ 
1 + — + 0{Hs 



1 + 2- 



N-2 c^M 



1 - - 
2 



5vr p^ 
1 iV-3 c^M 



ps 



+ 0{Hs) 



(2. 



Actually, these are first order results and — at least in the context of ferromagnetic 
systems — are well known [l^, ^ [ill]. We study these relations in Section 6 in order 
to fix the renormalization of the parameters. 



2.2 The free energy in "cubic" geometry {Lg = finite, T ^ 0, 
H, = 0,K = 0,l = 0(1)) 

We consider large volumes and low temperatures relative to typical non-magnon 
scales: 

-^<p,, T-^ps, (2.9) 
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while keeping the ratio / finite, 0(1). The free energy density reads 

T N -IT 



2 LI 



(2.10) 



Here qi and q2 are two non-universal constants which are not fixed by chiral pertur- 
bation theory. The constant qi is the infinite volume ground state energy density 
which is studied intensively for the model in eq. (|1 . 1| ) numerically 0, |20|, ^ |3^ 
and otherwise |^ |2^, (For a detailed discussion, see refs. |^8|, ^.) The 
constant q2 depends on the precise normalization of the partition function. The in- 
ternal energy density is already free of g2- The shape coefficients /3n(0 expressed 



through «„(/) by the relations [TT 



/5o(/) = «o(0 + ln(47r) - 7^ - - = ao(/) + 1.287142 



1 

An 



Oin{l) + 



n 



1,2, 



n{2n - 3) 

The functions ap{l) can be calculated easily using the following representation 



(2.11) 



ttp(0 = ap-3/2(0 + Q 



where 



with 



dir{l) 



dtf 
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IH 



S 



S{x) = ^e" 



(2.12) 



(2.13) 



(2.14) 



For / = 1 the relevant /3„ values are given in Table 1 of ref. |Tl|. In Appendix B we 
give approximate analytic forms, which are sufficiently precise for applications. 



The internal energy density e(T, Lg 
the form 



-T'^-§f ifiT,Ls)/T) can be written in 



e(T, Ls 



N -IT 

6 Z! 



(iV-2) 



he 

PsLgl 



d 



m)-ijim) +■■■ 



(2.15) 



In ref. [0 we quoted this result for = 3 with a slight errorQ, unfortunately: in 
the coefficient in the brackets the number " 1" was incorrectly given as " 2" . 



^Neither the text nor the subsequent equations are influenced by this error in 
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2.3 Thermodynamics at infinite volume (Lg = oo, T 7^ 0, 

H, = 0,hu = 0) 



Although the results in 2.2 were derived under the condition I =0(1), the ^ 00 
{I — * 0) limit is smooth^ and gives correctly the temperature dependence of the free 
energy and internal energy densities: 

/(T) = q,-(N-l)-^^T' + 0{T'), 



2n{hc) 

e(T) = q, + {N-l)-^^T' + 0{T'). (2.16) 



Here C(3) = 1.2020569. There is no ~ T contribution. Eq. ( |2.16D gives for the 
specific heat 

C{T) = (iV - 1)^^T2 + O(r^) . (2.17) 
nine)'' 

The ~ term in eq. ( p. 161) is just the free energy contribution of — 1 massless 



free bosons. This prediction has occurred earlier in [|T^ for the AF Heisenberg 
model (A^ = 3). Although our calculation goes up to order ~ T^, the corresponding 
coefficient turns out to be zero. There is an amazing mess in the literature on the 
~ term in the internal energy density of the AF Heisenberg model. We shall 
return to this problem in Section 3. 



2.4 Volume dependence of the ground state energy density 

(L, = finite, T = 0, Hs = 0, hu = 0) 

Eq. ( p. 101 ) was derived under the condition / =0(1) and as we shall discuss later, 
the T ^ (Z — > 00) formal limit of eq. ( p.lO| ) does not give the correct volume de- 
pendence of the ground state energy density. Very low temperatures require special 
care. We give two different derivations leading to 

N- 1 1.437745(^c) iN-l)iN-2) (hcf /l\ 

LI + 8 'JJA \^)' ^ ^ 

The ~ 1/L^ correction has been derived earlier for the AF Heisenberg model (A^ = 3) 
using current algebra like techniques |16| and effective Lagrangeans Jl^. In some 



works the ~ 1/L^ contribution to the ground state energy density was predicted 
to be zero, which is incorrect. The problem of the ground state energy has been 
discussed in the context of QCD in = 3 + 1 with SU(ra)xSU(n) symmetry Since 
SU(2)xSU(2)~0(4), eq. ( |2.18|) for A^ = 4 should agree with the corresponding result 
for n = 2 in ^j. After trivial adjustments corresponding to(i = 2 + l — > (i = 3 + l 
they agree, indeed. 



'see Appendix B 
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2.5 The energy of the lowest lying excitations and its vol- 
ume dependence {Lg = finite, T = 0, Hs = 0, hu = 0) 



It is known since a long time that the tower of lowest lying excitations in the broken 
phase is related to the slow precession of the order parameter and, in leading order, 



it is given by simple quantum mechanics over the symmetry group space [13, 21 



We have calculated the leading finite size correction also and obtained for the energy 
difference between the excited state with j = 1, 2, . . . and the ground stateQ 



E,iL,)~EoiL,)=jij + N-2) 



{hcl_ 



(2.19) 

For the AF Heisenberg model (A^ = 3) the leading order result has been discussed 
in 1116 . 



2.6 Weak staggered field in "cubic" geometry [Lg = finite, 
T ^ 0, 7^ 0, /i, = 0, / = 0(1)) 



We quote now the if^^dependent part of the free energy density derived under 
conditions of eq. ( ^.91 ) andQ 



Mo = 



^<0(l). 



(2.20) 



This result was derived earlier [|rT|, we restored the correct spin wave velocity de- 
pendence only: 



f{T,Ls,H,) = f{T,L, 



T 



ln(y.(n))+,..M%... 



(2.21) 



where the Hg independent part, f{T,Ls,Hs = 0) = f{T,Ls) is given by eq. ( p.lO| ). 
We also introduced the notations^ 



u = 7i'Uo , 



72 



A^- 1 
4 



2 PsLsl 



PsLgl , 



1(0^-2/52(0 , 



(2.22) 



^Eq. ( 2.19 ) refers to the energ y and not to the energy density: Eq{Ls) — L^eQ{Ls)- 
'^With the < sign in eq. ( 2.20 ) we want to indicate that Hs can be arbitrarily small. 



We have changed the notation pi, p2 of to 71, 72 to avoid confusion with the spin stiffness. 
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and Pp{l) was defined earlier in eqs. ( |2.11| - pTT^ . The function Yn{z) can be expressed 



in terms of tlie modified Bessel function Iu{z), 



I N-2 (z) . 



(2.23) 



Some of tlie properties of tlie function Yf^{z) are summarized in Appendix F of 



ref. 11 



Derivatives with respect to the staggered field Hs give the staggered magnetiza- 
tion and susceptibilities, 

2 



Ms{T,L,,H, 



d 



■f{T,Ls,Hs)=Ms{ii 



Y'Nju) 
Yn{u) 



+ 272M0 



he 



+ ... 



Xl{T, Ls, H,) = — = a 



Y"m(u) (Y'n{u 



Yn{u) 
Ms{T,Ls,H, 



Yn{u) 



55, 



+ 272 



he 

PsLJ 



+ 



(2.24) 



(2.25) 



In the Hg limit Ais{T, Lg, Hg) ^ as it should, and 



+ 



N -I ( he 



N' ''T 
2 



he 



(/?i(/)2 + (2iV-3)/?2(/))+... . (2.26) 



2.7 Volume dependence of the staggered susceptibility at 
small temperature (T <C hc/Lg, Hs = 0, hu = 0) 



As we mentioned before, low temperatures require special treatment. For T <^ he/ Lg 
the coefficients multiplying the expansion parameter he/psLg in eq. ( p.26[ ) become 
large. The dynamics is dominated by the low lying modes of the effective O(A^) 
rotator. The corresponding susceptibility is calculated in Section 8. We note here 
only that for temperatures 

TL'gPs , ^2.27) 



{hey 



< 1 



the susceptibility becomes independent of the temperature and it is proportional to 
the square of the spatial volume: 



Xs{T, Ls 



mIps 



N{N - 1) {hey 



1 + {2N- 3) 



he 3.900265 



PsLs 



4:71 



+ 




f2.28) 
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The basic reason of this strong volume dependence hes in the smallness of the 
excitation energies which are proportional to the inverse spatial volume, eq. ( p.l9|) . 



2.8 Uniform susceptibility {H^ = 0, hu = 0) 



For cubic geometry, T = 0(Tic/ Ls) and under the conditions in eq. 
for the uniform susceptibility at /i„ = 



2 Ps 
N {hcf ' 

+ — 



^ he 5 
1 + — —Jiil) 



3 psLsl 



we obtain 



(2.29) 



he 

^PsLsl 



1 ^ 



[N-2)m-^^Piiif-m) 



Here we introduced the shorthand notations 



(2.30) 



In the next -to-leading order a new shape coefficient ■?/'(/) enters, whose definition is 
given in Section 9. For / = 1 it can be expressed through 



-0.020529 



(2.31) 



for other values of / it requires a less trivial numerical evaluation. In Appendix B 
we give approximate analytic expressions for 



In the infinite volume limit eq. ( p.29| ) gives 



Xu{T) 



2 Ps 
N {hey 



l + (iV-2)-^ + (iV-2) 
z-nps 



T 



27rps 



+ ... 



(2.32) 



Here Xu{T) is the angular average of the transversal and longitudinal susceptibilities. 
This is refiected in the result Xu{T ^ 0) = (2/A^)p,/(^c)2 = (2/A^)x^. 

For T <^ he/Ls the result is again given by the effective O(A^) rotator. In 
particular for very small temperatures TL^ps/ (he)^ ^ 1 the susceptibility becomes 
exponentially small 



Xu{T, Ls 



2 1 



exp 



N - 1 {hef ' 
2 psL^sT 



(2.33) 
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3 Chiral perturbation theory vs. spin wave ex- 
pansion. Discussion on some problematic re- 
sults in the literature 



Chiral perturbation theory is an exact low energy expansion. The basic condition for 
the validity of this expansion is that the scales where the magnon physics is probed 
(like T or hc/Lg) should be much smaller than the typical non-magnon scales (like 
Ps). In eqs. (|2.1CI|J2.18| ), for example, hc/psLg enters as an expansion parameter. 
For large spin ps ~ JS"^, while he ~ JS, therefore the expansion parameter 
is he/ psLg ~ 1/5*. Since the spin wave expansion is an 1/5* expansion, for certain 
quantities the two expansions show analogies. The differences are more important, 
however. Consider eq. ( p. 16 ), for example. This result says that the internal energy 



of the AF Heisenberg model satisfies 

e(r)-e(0)^jCi3) 

which is an exact statement. There are no corrections to this result from higher 
order chiral perturbation theory or from other sources. Spin wave expansion gives 
only a l/S" expansion for the ~ coefficient. On the other hand, as it is discussed 
in the introduction, eq. ( p.lD is a symmetry relation only, it is satisfied for any S. 
Consequently, replacing the left and right hand sides of eq. (|3.1|) by the corresponding 
spin wave expansion one obtains an identity. This might serve as a good check on 
the results. 

As a further example consider the temperature dependence of the uniform sus- 
ceptibility in the large volume limit, eq. (|2.32 ). By replacing ps and he in it by their 



spin wave expansion form, which is available up to 0(1/5*^) [^, one obtains the 
spin wave expansion for the uniform susceptibility. This result is consistent with a 
direct spin wave expansion for Xu{T) up to the known order [Q. 

In the following we would like to discuss some problematic results in the AF 
Heisenberg model. 



3.1 The leading temperature dependence of the internal en- 
ergy 



The exact result, as we discussed before, is given by eq. ( p.l|) . Leading spin wave 
expansion predicts he = \/2Ja which gives on the right hand side of eq. ( |3.1| ) 
0.3826/( Ja)^. This is in agreement with the leading spin wave result for the left 
hand side Including next-to-leading corrections one obtains he = \/2Ji 



a 



'We consider = 3 and discuss the AF Heisenberg model in this section. 
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1.158. Inserting this number in eq. (|3.1| ), the ~ coefficient is predicted to be 
0.2853/( Ja)^. This is in agreement with Takahashi's modified^ spin wave theory, 



eq. (26) in ref. [|4 



These results are consistent. Schwinger boson mean field theory leads to 
the same equations as the modified spin wave theory in These equations were 
studied numerically in and (0.77±0.03)/(Ja)^ was predicted for the right hand 
side of eq. ( |3.1| ). This is almost a factor 3 larger than the correct result and is 

^ In 



due, presumably, to some numerical problems in solving the equations [26 



this incorrect result was confirmed in a Monte Carlo calculation. The reason might 
be that the simulation did not reach sufficiently low temperatures to be compared 
with a C(T) ~ behaviour. In |2^, the modified spin wave result for e(T) was 
incorrectly quoted as 2 ■ 0.2853/( Ja)^ • T^. This factor of 2 infiuenced some related 
conclusions in ref. [ES 



3.2 The 0(1/L^) correction to the ground state energy den- 
sity 



As we discussed in Section 2.5, the lowest excitations are described by a quantum 
mechanical rotator giving a ~ j(j + 1) ■ l/L^ contribution to the energy density. 
Although this term is zero for j = (ground state), there are non-zero ~ 1/L^ cor- 



rections to the ground state energy density from other fiuctuations (see eq. ( |2.18| )). 
This term is missing in earlier discussions and analyses [2^, |2^, For = 3 we 
can write 



lA377{hc) 



he 



5.7508p, 




(3.2) 



If we take ps/ J 
(4,. 



0.2, hc/Ja ^ 1.6 



29 



the correction is significant for Ls/a G 
12) where the best numerical results are available. A l/L^ fit might lead to 



a significantly distorted spin wave velocity. 



3.3 Uniform magnetic field at infinite volume and at low 
temperatures 

The uniform external field is coupled to a conserved quantity in the AF Heisenberg 
model. It enters the Hamilton operator like a chemical potential. In the effective 
Lagrangean this chemical potential shows up like the t-component of a constant 
imaginary non-Abelian gauge field. Here t is the third ("temperature") direction in 
the Lg X Lg X {h/T) euclidean box. 

-'^^ Takahashi's result corresponds to inserting he at S = 1/2 into the denominator in eq. 
rather than expanding the expression in 1/5. 



11 



Switch off now the magnetic field and take Lg = oo, T <^ ps- As it is observed 
before, new non-perturbative effects occur in this hmit, which can be related to the 
d = 2 classical non-linear sigma model |]15|, |T8|. The oo x (yo x (h/T) slab is reduced 
to an oo X cxD two dimensional model. The intuitive reason is that the correlation 
length along the spatial dimensions is exponentially large in 1/T, so the width of 
the slab (h/T) becomes very small compared to the correlation length. In the d = 2 
non-linear sigma model it is an excellent technical tool to probe the system under a 
chemical potential corresponding to a constant gauge field along one of the two spa- 



tial directions |T8[ . Among others, exact results can be obtained for the chemical 
potential dependence of the free energy. This chemical potential has, however, no 
direct relation to the uniform magnetic field in the AF Heisenberg model, which 
corresponds to a constant gauge field along the "temperature" direction. This fact 
modifies the conclusions of a recent paper significantly. 



4 The effective Lagrangean. Discussion on the 
cut— off effects 



In this section we construct the effective Lagrangean which serves to derive the free 
energy density up to next-to-leading order. Since for quantum field theories or 
critical statistical systems this problem has been discussed in detail before [^, 0, |^, 



11], we can be brief on those points. 



The AF quantum Heisenberg model in eq. ( |1.1| ) has two parameters: the coupling 
J which sets the scale and the value of the spin 5* which takes discrete values. The 
smaller is S, the closer is the system to the boundary between the broken and 
symmetric phases. It is rigorously established that the ground state is ordered for 
S" > 1 [^. Approximate analytic and numerical results suggest that = 1/2 is in 
the broken phase also For every S we have two massless excitations (magnons) 



with a relativistic dispersion relation. This is similar to the case of a continuum 
quantum field theory, or a classical ferromagnet in the critical region with a broken 
0(3) symmetry. On the other hand, in the later cases the non-Goldstone boson 
scales (e.g. ps) are also much smaller than the cut-off. This is not so in the model 
in eq. ( |1.1| ). Take, for example the length scale corresponding to the spin stiffness, 
= he/ Ps- For an order of magnitude estimate we can take he = 1.6 ■ Ja, ps = 0.2 ■ J 
|]28| , |29[| giving ^p/a = 8, as opposed to a critical system where this ratio goes to 
infinity. In general, this cut-off dependence modifies the effective prescription. 

In the method of effective Lagrangeans one writes down the most general local 
Lagrangean Lg// (in terms of the smooth field corresponding to the order param- 
eter) which respects the symmetries of the underlying model. The natural — al- 
though strictly unproved — expectation is that the predictions from Le// reflect the 
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symmetries only and correspond to any underlying model having the same kind of 
Goldstone bosons (magnons) and symmetries. This method becomes practically in- 
teresting through the observation that Goldstone bosons (magnons) interact weakly 
at low energies and a systematic perturbative expansion can be advised. In any 
given order only a finite number of unknown low energy constants enter. 

Forget about the cut-off effects first and consider the problem of the free energy 
density /(T, Lg) for small T and under the condition / =0(1) where / is defined 
in eq. ( p.3| ). Let us take "he = 1" in the following discussion, the missing he factors 
can be inserted at the end using dimensional analysis. The effective Lagrangean 
is constructed as a general non-linear cr-model in terms of the field variables S, 
S2 = 1: 

S(x) = ((l-7r2(x))'^',7f(x)) , (4.1) 
TT (x) = (n^ix), 7r^(x)) , = {xq = t, Xi, X2) . 

The low energy excitations carry momenta p ~ Lj^ ~ T. Then, in the effective 
Lagrangean, every derivative 9^ is counted as ~ p. The terms in the effective model 
should comply with the symmetries (most notably with the 0(3) symmetry) of the 
underlying theory in eq. ( |1 . 1| ) . The different terms in the effective Lagrangean are 
multiplied with unknown couplings. The leading term in the effective action contains 
the least number of derivatives[^, 

A = f'^ dt dPx \psd^S{x)d^S{x) (4.2) 

t/ «/ ^ 

, /"^^ ,2 1 fo - (7fa7r)(7f9.7f)l 

dt / d^x -ps d^nd^n + ^ ^ ' , /i = 0, 1, 2 

JO 2 1 — TT^ 



In eq. (t4.2|), we denoted the coupling by ps using the well known result that for 
T = 0, = 00 this coupling is the spin stiffness (hehcity modulus, square of the 
Goldstone boson decay constant) and the fact that the couplings of the effective 
Lagrangean are independent of T and Lg 0.0 The field vr should be counted as 
~ p^''^, since fiuctuations of this size have a Boltzmann factor of 0(1). Indeed, the 



leading term in the Lagrange density in eq. (4^) is 

1 



^Psd^TTdf^n p^ , (4.3) 



which is integrated over a region L^/T ~ 1/p^. The leading term in eq. ([4.3|) gives 
an O(p^) contribution to the free energy density. The terms obtained by expanding 
the denominator in the second term in eq. ([4.2| ) give contributions of O(p^), O(p^), 



repeated indices are summed 
^^This is true if periodic boundary conditions are used which we shall assume in the following. 
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. . . . Of course, there are other terms in the general effective Lagrangean which can 
also contribute in higher order. The terms with four derivatives have the form 



\9^l^ {d,d,Sd,d,S) + \gf^ {d,^d,^f + \gf^ (9,S9.S) (a^Sa.S) . (4.4) 

By a change of integration variable in the path integral the first term in eq. ( [4.4| ) 
can be transformed away ITlI], so we can take g^l^ = 0. As it is easy to see, the other 



terms in eq. ([4.4|) are at least 0{p^). Terms in the effective Lagrangean containing 
more then four derivatives are even more suppressed. We can conclude therefore that 
if the cut-off effects can be neglected, the effective action in eq. ( [4.2|) determines 
the O(p^), O(p^) and O(p^) parts in the free energy density. Only two unknown 
parameters, ps and — after restoring dimensions — c, enter up to this order. 

Let us now discuss the cut-off effects. The effective prescription is influenced by 
a finite cut-off in two ways. First, the finite cut-off might break some of the sym- 
metries (for example, space rotation symmetry) which the underlying model would 
have in the continuum limit. This results in new terms in the effective Lagrangean 
which at infinite cut-off would be excluded. Second, the Fourier integrals (sums) of 
the chiral perturbation theory will depend on a finite cut-off. The divergent cut-off 
dependence is absorbed by the couplings through the usual renormalization proce- 
dure. There remain contributions, however, which contain inverse cut-off powers 
and depend on T and Lg. 

We want to argue now that the cut-off effects enter the free energy density on 
the 0{p^) level first. At finite cut-off spatial rotation symmetry is lost, only the 90° 
discrete symmetry of the square lattice remains. We need at least four derivatives 
to observe the difference. Since every term contains at least two tt fields, the new 
term is 0{p^) or higher. Actually, there is a new contribution at O(p^) already. The 
reduced symmetry allows the new four-derivative term 

Igt^ E (dASdAS) , (4.5) 

^ i=l,2 

which is invariant under 90° spatial rotation but not rotation symmetric. The cor- 
responding leading term 

Igt^ E (dAimdi^r) , (4.6) 

^ 1=1,2 

is O(p^). The leading contribution to the free energy density coming from this term 
is 

^Ly^ilA (47) 

where the denominator is the magnon propagator in leading order. It is easy to 
show that eq. ([4.7|) has a non-zero, 0{p^) temperature and volume dependence. 
Consider now the second source of cut-off dependence generated by the momentum 
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sums in chiral perturbation theory. The leading 0{p^) contribution is a 1-loop graph 
corresponding to a momentum sum.0 The leading cut-off effect in bosonic sums 
and integrals is ~ (l/A'^"*)^. We expect therefore ■ p^/(A'^"*)^ =0(p^) cut-off 
dependent corrections. 

These considerations show that the effective action in eq. ( [4.2| ), which after restor- 
ing the dimensions has the form 

f^/T , fLs 1 

/ dt d x-ps 
JO Jo Z 

is sufficient to derive the leading (0(p^)) and next-to-leading (0(p^)) contributions 
to the free energy density. These contributions depend therefore on two parameters 
Ps and c only and are cut-off independent. Eq. ( [4.8[ ) can be generalized trivially to 
O(A^) symmetry by defining S as an iV-vector and vf as an (A^— l)-vector. In the 
cut-off dependence has been investigated in the large- limit of the classical O(A^) 
ferromagnetic Heisenberg model. The exact solution showed a cut-off dependence 
which is consistent with our discussion above. 



=1,2 



(4. 



5 The free energy density /(T, Ls) up to next— to- 
leading order ("cubic" geometry) 

We take " ^ = c = 1" and restore the dimensions at the end. It is useful to introduce 
the notations 

Lt = -^, V = Ll-Lu L = V^/\ x = {xo = t,xi,X2). (5.1) 
In the effective prescription the partition function is given by 

Z = ll [ dS{x)6 (S\x) - l) ■ e-^(^) , (5.2) 

where the action A is given in eq. (f4.2|) . As we discussed before, the field 7f(x) is 
small, 0(]9^/^) which is the basis of the perturbative expansion in eq. (|5.2|) . The 
euclidean box X Lg X Lg is large compared to the scale of non-magnon excitations 
and we shall assume that its shape is essentially cubic 

>-, Li>-, fi = ^ = 0{l). (5.3) 

Ps Ps J^s 

^•^The summand is essentially the logarithm of the leading magnon propagator. See Section 5. 
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d xS(x) 



(5.4) 



Perturbation theory in a box of this kind has a technical problem due to the existence 
of zero modes. The physical reason is that the dominant S-configurations have a 
net (staggered) magnetization, i.e. a non-zero value for 

1 

V 

which freely rotates around in the O(A^) group space. One has to separate these 
modes before starting a perturbative expansion. Since the technique is standard we 
quote the result only |lT|. Eq. ( |5.2| ) can be written in the form 

z = n I / _ ^,^^^^1/2 ) n ^ / d'^^'i^) ) X 



i=l 



exp|-^(S) + (iV- l)ln 
Going over to Fourier-space, we have 

iV-l / 1 ^ 

J3 



— / d^x (1 - Ti'^ix 



1/2 



n / 



N-l 



n 5 / d'xA^)) - n (/ rf^i) n ^ (^u) ■ 



■ \nV 



(5.5) 



(5.6) 



Eliminating the zero modes therefore has the following consequences: in the Fourier 
sums the k = (0, 0, 0) mode should be left out, the action receives an additional, 
field-dependent contribution (the second term in the exponent in eq. ( p.5| )) and a 
field independent, but size dependent constant ^(iV — 1) liaV occurs in the exponent. 

One follows the rules of standard perturbation theory except that in calculating 
Feynman graphs in the momentum sums the zero mode does not enter. According 
to eqs. (|5.5|J5.6|) in leading and next-to-leading order the Lagrangean has the form 



C 

Co 



^0 + ^1, 



N-l 
2V 



\nV, 
N-l 



(5.7) 



-TT 



2' ' ' ' ' 2V 

Here £1 is 0(p'^) and can be treated as a perturbation. As discussed in Section 4, 
there are no T- and L^-dependent cut-off effects up to the order of our calculation, 
0(p^). Therefore, the results are universal, independent of the regularization used 
in the effective theory. A conventional possibility is dimensional regularization |^ . 
In this case the factor (1 — 7r^(a;))~^/^ in the measure in eq. (|5.5|) can be suppressed 
].Q One obtains for the free energy density 

/ = 
fo 



/0 + /1 
A^- 1 



2V 
fi = (^i)o 



A^- 1 



2V 



InV 



(5.8) 



^^Dimensional regularization is very convenient, but certain aspects of it are rather formal. We 
have checked that the final results of this section remain unchanged if lattice regularization is used. 
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where J2' denotes a summation over the discrete momenta fc^ = (fco, ki, of the 
periodic box LfXLgX Lg with the fc^ = (0, 0, 0) mode left out, while A;^ = kQ + kf + k"^. 
The notation ( )o refers to the expectation value with the leading Cq in the exponent: 



1 



1 



(£1)0 = —{N - l)G(O) OpyCix -y) + - l)^G(O) , (5.9) 

where 

_ 1 . pikx 

Gi^) = yj:,^- (5.10) 
In calculating the momentum sums in eqs. (|5.8|j5.9|) one separates first the V = 00 



part, which is in general divergent. This part is absorbed through the renormaliza- 
tion of the parameters of the effective Lagrangean. In dimensional regularization 
this renormalization is trivial, since the corresponding integrals are defined to be 
zero.|3 The remaining finite, T- and Lg-dependent part can be written in a form 
which allows an easy numerical evaluation for any value of the shape parameter 
= Lt/Lg. Since these steps were repeatedly discussed in the literature [0, ||, [TT| , 
we quote the final result only. Eqs. ( |578| , |5^ ) can be written as 



Jo = gi + 77 - 



V 2V 



where qi and q2 are two non-universal constants, / and the functions (3n{l) were 
defined in eq. ( |2.3| ) and eqs. ( p.ll| - |2TT^) , respectively. By restoring the dimensions. 



eq. (|5.11|) leads to the result quoted in eq. (|2.1CI| ). 



These results have been derived under the assumption that I = {hc/TLsY^^ is 
0(1). What happens if we take the T — >■ 0, or — > 00 limit? 

The T ^ {L ( — i> 00, or / — i> 00) limit corresponds to an infinitely elongated 
cylinder geometry. As it is well known |TB| , PT| , ^ in this geometry there are additional 
quasi-zero modes corresponding to the slow rotation of the (staggered) magnetiza- 
tion of the Lg X Ls planes as we move along the t-direction. These quasi-zero modes 
should be treated separately. The expansion used to derive eq. ( [5.11| ) breaks down 



in the I —* 00 limit and does not give the correct ground state energy. For later 
reference we quote the result of this limit, nevertheless: 

N -Ihc 

eq. (5.11)|,_ : fo —j^ ■ 1.437745, 

, , {N - 1){N - 2) jhcY 



^^In Sections 6 and 8, where lattice regularization will be used, we carry through this renormal- 
ization process explicitly. 
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There is an other way to see the problem with the T ^ hmit. In this hmit we 
have effectively a. d = 1 non-linear a-model. In d = 1, the Lf ^ oo limit can not be 
interchanged with the weak coupling expansion |^. Since this problem is specific 
to d = 1, we do not expect similar problems in the — »• cxd (/ — > 0) limit, where 
we obtain a two-dimensional slab. Although there are new non-perturbative effects 
in this limit which generate a finite mass m to the magnons |jl5|, |1^, this mass is 
exponentially small in 1/T. The effect of such a small mass is negligible in the free 
energy density assuming that there are no infrared divergences in the m — limit, 
i.e. assuming that the perturbative expansion for / in d = 2 is infrared finite. This 
is actually the case [Q. Eq. ( |5.11D gives then 

/(T) = gi - (iV - l)^T' + O (r^) . (5.13) 

There is no ~ contribution, since /i|;^o ~ ^- This is the result (after restoring 
dimensions) announced in eq. (|2.16|) . 



6 Volume dependence of the ground state energy 
and of the lowest lying excitations 

The special properties of the cylinder geometry corresponding to this problem have 
been observed and discussed before Jl^, ^ |]. The next-to-leading volume correc- 
tion to the ground state energy has been calculated in d = 4 for the symmetry group 
SU(n) xSU(n) relevant for the strong interactions The next-to-leading ground 
state energy for O(A^) and the volume correction to the lowest excitation energies 
to be discussed here are new results. 

Since this problem is a somewhat non-trivial application of chiral perturbation 
theory, we use this occasion to discuss some technical issues: we consider a "phys- 
ical" regularization (lattice, as opposed to dimensional regularization of Section 5) 
and discuss the renormalization process explicitly, especially the problem of power 
divergences. We give also some details on the way collective coordinates are intro- 
duced in treating the zero and quasi-zero modes. 

We shall consider a (i-dimensional euclidean cylinder Lg~^ x Lt in the Lt ^ oo 
limit. The AF Heisenberg model corresponds to ci = 3, = 3 in the following 
expressions. We shall use a (hyper) cubic lattice to regularize the effective theory. 
We put the lattice unit equal to 1 (a = 1) and denote the lattice points by a; = {t, x), 
where t = 1, 2, . . . , Lj, etc. The Fourier transform is defined as 

7i{x) = ^Y^^^'^m , V = Li-'-L,, (6.1) 

and we introduce the notation A^: 

^t^9{x) = g{x + fi) - g{x) , (6.2) 
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where fi is the unit vector in the /^-th direction. 



6.1 Introducing collective coordinates 

Until Lt =0{Ls), it is sufficient to take special care of the freely moving total 
magnetization. For large Lt/Lg, however, the magnetization of distant time slices in 
the cylinder can differ significantly. We have to treat these slowly moving variables 
(they are the k = {ko, k = 0) modes) non-perturbatively. 

The action in eq. ( |5.2|) has the regularized form 

^(S) = ^P°E^MS(x)A,S(a:), (6.3) 

X,fl 

where we denoted the coupling by anticipating a non-trivial renormalization. A 
possibility to introduce appropriate collective coordinates is to insert into the path 
integral 

r I- ^-1 

n 

/ 71=0 

where V, = L'i~^. We write 



(6.4) 



m(t) = m{t)e{t) , e{t'f = 1 , 

dm{t) = m{t)^-^dm{t)de{t) , (6.5) 

f , 27r^/2 
/ de = K = — ; — — . 
J T{N/2) 

Introduce the O(A^) rotation Q(t) 

e{t) = Q{t)n, n= (1,0,...,0), 

neO{N), Jdn=l. (6.6) 

Clearly, the vector e does not completely fix the matrix Q. On the other hand, 
integrating over Q rather than e gives an e-independent constant factor only. The 
partition function has the form 

Z = H [ dS{x)6[s\x) - l) n / dm{t)m{tf-^ x 

X •' t '' 

K J dn{t)5^^^ (m{t)n{t)n ^) j exp {-A (S)} . (6.7) 

Replace S by the new integration variable R: 

S{t,x) = n{t)R{t,x) . (6.8) 
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Since the measure is O(A^) invariant, we get 

Z = Y[ f dR{x)6 (r\x) - l) n / dm{t)m{t)^-^ x 

K J dn{t)6^^^ (m{t)n -yJZ ^(t, x) j exp {-A (fiR)} . (6.9) 

Write 

R^ = Tr\ i = l,...,N-l, R^=(l-TT^y^\ (6.10) 
and integrate over the 5-functions 

Z = n / d7r{x)l[ f dQit)l[6(^Y.^\t,x)] ■ exp {-A^iQ,n)} , 

X t i=l ^ I 

A(^^,vf) =^(f]R) + ^5]ln(l-7r2(x)) - (6.11) 

^ X 



t \ ^ X J 



Going over to Fourier variables, eq. ( 6.11 ) leads to 



\ k i=l J \ t '' / \kO i=l 



X 



N -1 

exp \-Ai{Q,n) - Lt^— In VA . (6.12) 



6.2 Integrating over the collective variables 

The action ^(R) can be written as 

^(R) = J2 + A) - = dVp'l - p° E + A) • R(^) • (6.13) 

x,^ x,ii 

Since fl is independent of x, we have 

d-l 

AinR) = dVp^^ - p° ^ ^ R(x + i)R{x) -pIY. ^(^ + + l,x)- n{t)R{t, x) . 

a; i=l t,x 

(6.14) 

Consider the sum in the last term of eq. ( |6.14]) 



^ [fi(l)^fi(2)R(2,f) ■R(l,f) +^](2)^fi(3)R(3,f) ■R(2,f) + ...] . (6.15) 

X 

As eq. (|6.6|) shows, the direction e{t) of the magnetization of a time slice is not 
influenced by the part of VL{t) which leaves the vector n invariant. These are spurious 
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degrees of freedom in the O(A^) matrix Q. We shall get rid of these spurious variables 
by changing integration variables in eq. (|6.12|) appropriately. Replace the integration 
variable n{2) in eq. (|6l2|) by 



V{2) = Q{lfQ{2), 
and introduce the rotated vectors 

R(2,f) = V{2)R{2,x) 
Here V{2) is an O(A^) element with the structure 
/ 1 ■■■ \ 



(6.16) 



(6.17) 



y(2) 



■• 



i,k = 1,2, 



AT- 1 



(6.18) 



The 0(A^ — 1) rotation Vik{2) serves to exclude the spurious degrees of freedom in 
V{2) and will be specified further later. Eq. ( |6.17D is an 0(A^ — 1) rotation on the 
integration variables 7r*(2,a;) 



7r^(2,f) = V,,(2)7r'=(2,f). 



(6.19) 



We can introduce 7r*(2, x) as new integration variables in the t = 2 plane in eq. ( |6.12| ). 
The measure is invariant and the 6 (^7r^{ko,k = 0)^ constraint is as well since the 
rotation is independent of x. The terms in the exponent in eq. ( p.l2| ) are also 
invariant, except that eq. ( |6.15|) goes over to 

r(2)f(2)^R(2,f) ■R(l,f) + (^(l)y(2)y(2)^)^^](3)R(3,f) ■R(2,f)... . 

(6.20) 

Now we can repeat this procedure with the second term in eq. ( |6.20| ). We replace 
the integration variable Q{3) by 

V{3) = (^(1)^(2)^(2)^)"^ n{3) , (6.21) 

and transform the variables in the t = 3 plane by the x-independent rotation 

R(3,f) = t>(3)R(3,f), (6.22) 

where V^(3) has the structure of eq. ( |6.18| ). Repeating these steps in the t = 4, . . . 
planes, eq. (|6.2CI| ) goes over 

J2 ['^(2)V^(2)^R(2, x) ■ R(l, x) + ■\/(3)V(3)^R(3, x) ■ R(2, x) + . . .] . (6.23) 
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Boundary effects can be neglected for Lt ^ oo. Introducing the notation 

V(t - 1) = V{t)V{tf (6.24) 
and suppressing the "hat" on R, we get 

z=(iij dnk)^ (u J dv{t)^ mi[s{nk',k = 0) 



exp-|^(y,R) + i^ln(l-7f2( 



x)) - (6.25) 



(iV-l)^ln(^l^i?0(t,f)^ -Liln/t-L,^^lnK| , 

where 

X i=l t,x 

= ^P°E^MR^MR-P°EQWI^(^+1'^)-I^(i'^)' (6-26) 

t,x 

with the notation 

Q(t) = V{t) - 1 . (6.27) 
We write the exponent in eq. ( |6.25| ) as 

{,..} = A^-) + A^-y) -Lt\nK-Lt^^^^\nVs, (6.28) 

where ^'^^ is the vr-dependent part, while ^'^'^^ is the V^-dependent part including 
the TT-V interactions. Since we have replaced the dangerous components 7f*(/co, k = 
0) by collective variables, we can consider the remaining vf variables as small, and 
expand. Eqs. ( |6.25| , |6.26D give for the vr-dependent part 



4^^) = lp°5:A,7fA,7f, (6.29) 

X,fl 



4^^ = ^p°E 



^(7fA^7f)(7rA^7f) + (7fA^7f)(A^7fA^7r) + 

X,fl 

\A^7fA^7f)(A^7fA^7r) 



4 



N -l^^o 1 



OT/ ^ 



2K V 2 



The last term in Ai^^ comes from the measure, it plays the role of a counterterm in 
the renormalization as we shall see. In deriving eq. ( |6.29| ) we used the identity 

A^ {f{x)g{x)) = A^f{x) ■ g{x) + f{x) ■ A,g{x) + AJ{x) ■ A,g{x) , (6.30) 
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(no summation over /i). The \^-dependent part is generated by the last term in 
eq. (|6.26|) . Expanding in vr, we get for this term 

- E {Qooit)p% - Qooit)A{t) + Q^kit)B,kit) + Q^kit)Qkit) + 0(7r7r7r)} , (6.31) 

t 

where 

Mt) = p°Ek(^'^) + ^^o(7f'(t,x))j , (6.32) 



Bikif) 

Cikit) 



7fi{t,x)7Ck{t,x) + -Ao7fi{t,x)7ik{t,x) + -TTj (t , f) AoTTfc (t, f) 



-p° E [^07fj(t, X)7ffe(t, X) - 7fi(t, x)Ao7[k{t, X)] . 



As eqs. ( p.25| , |63ID show, in every time shce t we have to perform the group integral 
J dVexp {p°Lf-i (V^oo - l) - QooA + QikB.k + QikQk} , (6.33) 

where V and Q are related to V through eqs. ( |6.24| , |6!27D . Due to the form of V, 
eq. (|6.18|) , we have 

Voo = Voo , Qoo = Voo-l. (6.34) 
For p^L'^~^ large, the typical V matrices contributing to the integral in eq. (|6.33|) 

1 \ 



are 



(6.35) 



where e is small. Then V^q =C)(^)- We shall now specify the matrix V in eq. (|6.18 



in terms of the matrix elements Vik in such a way that the integrand in eq. ( H.33 ) 
becomes free of the spurious degrees of freedom discussed before. We shall take V^fc 
essentially equal to Vik, plus corrections needed to make Vik an element of 0(A^ — 1): 



For small e we have 



Vik = OiiVik + E Aj^jfc , no sum over 
i 

f^ij /3ji , /3ii . 



= 1 + ^ {Vof + 0{e^) 



(6.36) 



(6.37) 



f3i, = ^VoV,o + 0{e'), t^j. 
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Then V in eq. ( |6.24| ) has the form 



V = VV 



iO 



1 - 



(6.38) 



Only Voo and Vio, with Vqq + Z^i ~ ^ enter in eq. ( |6.38|) , representing the motion 
of e in eq. (|6.6| ) without additional spurious degrees of freedom. Eq. ( |6.38 ) shows 
that Qoo ~ Qik ~0(l/p°Lf^^). This is just the order we want to go. Expanding 
eq. (|6.33|) in terms of Qoo and Qik, we get 



dVexp{p^^Li-'Qoo} 



1 — QoqA + ^ QikB: 



i,k 



i,k 



where 



Q 



00 



K)0 ~ 1 ) Qik 



iOVkO ■ 



(6.39) 



(6.40) 



In eq. ( |6.39|) we can replace A and Bik by their expectation value calculated with 
the leading part of A^'^'> (the first term in eq. ( |6.29| )): 



where 



A^{A) = {N -l)VsD*{Q), 

B^k ^ 5,k{B) = 5,kVs {D*{{)) + AoD*(0)) 



In eq. ( |6.42| ) a constrained summation over the momenta occurs, 



Using 



N-l 



i=l 



« k^{kO,k=0) 
1 N-l 1 

4 E = -l{l- Vol) = Qoo + 0(e^ 
i=i ^ 



we can write eq. ( |6.39 ) in the form 

1 dV{t) exp {ee//Qoo} , e^// = p°Lfi - {A) + (5) . 



(6.41) 



(6.42) 



(6.43) 



(6.44) 



(6.45) 



At this stage it is clear that eq. ( |6.45| ) represents the motion of an O(A^) rotator 
with inertia 6e//- Indeed, we have 



Qoo(t) = K,o(t) - 1 = {^{t - ^fm))^^ - 1 = fi(t _ l)n . fi(t)n - 1 



(6.46) 
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and the integral in eq. (|6.45| ) can be written as 

^ / deit) exp [-l^eff (e(t) - e(t - 1))'} . (6.47) 

Here we used that dV{t) = dQ{t) = de{t)/K. 

The group integral in eq. ( |6.45|) can be easily performed. For large p^V^ we get 

where 



^ ;N\ 1 N-i 



2 / 2^ 

The corresponding part of the free energy density reads 



/ 



(7r,V) 



A^- 1 

N -1 



(iV- l)(iV-3) 



[{N -2)D*{0) - AqD* 



(6.49) 



(6.50) 



6.3 Contribution from the perturbative fast modes 



The leading contribution is given by the first term in eq. ( |6.29| ) which leads to 



Jo 



(n) 



N -1. 
2V 



where 



d{k) = A J2sm^'^ 



(6.51) 



(6.52) 



and Yl* is defined in eq. ( |6.43|) . We can write 



J2*^^{d{k))= J2 ln(rf(A;)) - ^ ln(rf(A;°,A^ = 0)) . (6.53) 

k kyt{0,...,0) kOjtO 

The sums on the right hand side of eq. ( |6.53D are standard in chiral perturbation 
theory and can be evaluated easily. In the ^ oo limit one obtain^ 



N -1 
2Li 



a%']{l) + 2-'- 



2Lf-i 27r 



(6.54) 



where are the "shape coefficients" defined earlier for d = 3 (eqs. ( p.l2| - |2rT^ )). 
In eq. (|6.54|) we need the shape coefficients of a (d — l) dimensional cube only which 



16 



We suppress the non-universal terms ^const. or ^^const./y in /. 
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are listed in Table 1 in ref. [^. The next-to-leading correction is given by the 
expectation value of in eq. ( |6.29| ). Using the relations (Appendix A) 



-2A^Z}*(0) (no sum over /i) 



(6.55) 



one obtains 



/i 



(^) 



Pi 



^A,D*{0)A,D*{0) + ^^^^ ^^ D*(0) 



(6.56) 



6.4 The ground state energy 



Collecting the size dependent constants (the last two terms in eq. (|6.28|) ), the con- 
tributions from the collective coordinates (eq. ( |6.50| )) and from the fast modes 
(eq. ( |6.56|) ) we get a surprisingly simple result for the ground state density 



eo(^s 



lim /(T,L.) 



N 



2Li 
N -I 



a 



{d-i). 



2p° 



V.(l) + 2-- 
A^D*(0)A^D*(0) 



(6.57) 



Ao/^*(0) 



+ 



(A^- l)(A^-3) 



In obtaining eq. ( |6.57|) we used the definitions of n and u in eq. ( |6.5|) and (|6.49| ), re- 
spectively. In order to proceed further we need the explicit form of the lattice Green 
functions in eq. ( |6.57|) . As it is discussed in Appendix A, all the dangerous cut-off 
dependence cancels in eq. ( |6.57| ), and up to a non-universal volume independent 
constant we can write 



A^Z}*(0)A^Z}*(0) - -^Ao/^*(0) 



(6.58) 



The spin-stiffness occurs in eq. (|6.57|) first in the correction only, we can therefore 
replace p° by the renormalized, physical ps- (The renormalization properties will be 
discussed in 6.7). We get then {Vs = L'^~^) 



eo(^. 



N 



eo[oo) 



2Li 



a 



{d-i) 

-1/2 



(1) 



2' 

^-1 



+ 



{N -1){N -2) 



(6.59) 



The shape coefficients a"^^l^{l) are listed in Table 1 in ref ITT 



The result quoted in 

eq. ( |2.18| ) is obtained after restoring the dimensions and setting c? = 3. For = 4, 
(i = 4 eq. (|6.59|) agrees with the result in ref. p obtained for SU(2) xSU(2). 
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6.5 An alternative, intuitive method 



Since the previous calculation is rather involved, it is reassuring that the result in 
eq. ( |6.59D can be reproduced in a simple intuitive way. The following considerations 
offer additional insight into the physics of the problem also. 

In Section 5 we studied the free energy density /(T, Lg) under conditions T -C 
Ps, l/Ls <^ Ps, TLg =fixed constant. Very low temperature corresponds chosing 
this fixed constant very small, but even in this limit the temperature T remains 
~ l/Lg. Standard Goldstone boson (magnon) excitations carry momenta ~ l/L^ 
and they freeze out at these low temperatures. The rotator excitations, however, 
carry energies ~ and for those excitations the temperature is high even in this 
limit. They are fully excited and included in the free energy. This is the reason why 
eq. (|5.12| ) disagrees with the correct ground state energy density eq. ( |2.18|) . Rather, 
eq. (|5.12|) should agree with the free energy density coming from the O(A^) rotator 



H 



1 



26 



eff 



(6.60) 



in the high temperature limit, where Oe// ~ PsLi and Eq{Ls) is the ground state 
energy we are looking for. The operator is, in coordinate representation, the 
Laplace operator on the A^ — 1 dimensional unit sphere. This matching condition can 
be used to calculate the ground state energy Eo{Ls). We need the high temperature 
expansion of the partition function 



7rot 



Trexp — 



IT 
T 



exp 



T 



Tr exp 



2ee//T 



(6.61; 



An elegant way to obtain the result is to use Seeley's expansion giving 



Tr exp ( — L t 



^^-(JV-l)/2 



^ ^ (iV-l)(iV-2) ^ 
6 



(A^-l)(iV-2)(5iV2-17iV+18) 2 
^ 360 ^ 



(6.62) 



with the normalization factor A = ^^"^ /V{N /2)^ which will not be important 

for us. The free energy density predicted by the rotator problem has the form 

Matching eq. ( |6.63|) with eq. ( |5.12|) we find the next-to-leading correction to eo{Ls) 



{N-l){N-2) {N-l){N-2) {N-l){N-2) 



+ 



12p,K' 24p,K' 
which is the result obtained before, eq. ( |6.59|) . 



(6.64) 
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6.6 The lowest lying excitations and their volume depen- 
dence 



The collective slow modes correspond to an O(A^) rotator with inertia G ~ -Z^f"^- 
The corresponding excitations are ~ l/L^^^, lying much below the excitations of the 
fast TT-modes which are ~ ^/Lg. Using eqs. ( |6.40 , 5!^ , 6.44 ) we can write eq. (|6.39|) 
in the form 

j d\/exp{(pOLfi - {A) + {B)) (Voo - 1)} • (6.65) 
Eq. (|6.65|) shows that the 7r-modes modify the inertia of the O(A^) rotator to 



Td-l 

s 



{A) + (B) = V, 



(N -2)D*{0) + AoD*{0) 



(6.66) 



As we shall discuss in 6.7, eq. (|6.82|) , the relation between the bare and renormalized 
spin stiffness in our order reads 



p° = + (iV - 2)Doo(0) + 



1 

2d 



(6.67) 



where Doo{x) is the infinite volume lattice propagator. Using the explicit form of 
the Green functions (Appendix A) we get 



Qeff = PsL 



d-l 



1 + 



-a 



(d-i) 



The excitations of the O(A^) rotator are given by 



(l) + 2 



d 



d-2 



(6.68) 



Ej — Eq 



26 



eff 



(6.69) 



which, at (i = 3, leads to the result quoted in eq. (|2.19|) . 



6.7 Renormalizations. Relations at infinite volume and zero 
temperature 

The physical parameters ps, c and M.s are defined at infinite spatial volume and zero 
temperature. In the following we study briefly their renormalization with lattice 
regularization in ci = 3 and, as a side product, we rederive the results quoted in 



eq. (13) 



Switch on an external symmetry breaking staggered field = {Hs,0, . . . ,0). 



The corresponding leading term in the effective action reads ITT 



5:A4°H,S(x), (6.70) 
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where is the bare staggered magnetization. We consider a d{= 3) dimensional 
infinitely large hypercubic lattice. There are no special zero modes in this case. 
Expanding the action (eq. (|6.3|) and eq. (|6.70|) ) and the measure in the vf fields one 
obtains 



where 
Co = 

Ci = 



Z = Y[J rf7f(x)exp|-^(/:o + A + 



(6.71) 



[7fA^7r)(7fA^7f) + (7fA^7f)(A^7rA^7r) + -(A^vf A^7f)(A^7rA^7r) 



with the notation 



-71 

2 



Let us calculate first the (staggered) magnetization: 



M"" 

VZdH. 



(6.72) 
(6.73) 



^^L^'-HO)). (6.74) 



Eq. ( |6.74| ) and the subsequent equations of this section are understood as relations 
in the first non-trivial order. In eq. ( |6.74| ) D^°\x) is the infinite volume massive 
lattice propagator 



X] 



Akx 



(6.75) 



{2TiYd{k) + M^ ' 

where d{k) is defined in eq. (|6.52|) . At Hs = 0, goes over to the spontaneous 
magnetization Ms which determines the relation between the bare and renormalized 
magnetizations 



Ms = M'! 1 



2ps 



(6.76) 



where we replaced p° by ps since the difference is of higher order. Eqs. (|6.74| , pr7BD 
and the relation 



D(^^«)(0) = dW(o)_ 



Mo 



lead to the result 



M 



H 



Ms 1 + 



47r 

M 



{d = 3) 



57r 



Ps 



(6.77) 



(6.78) 



where we replaced Mq by M = AigHs/ Ps- After restoring the dimensions, eq. (|6.78 



gives the result quoted in eq. (|2.8|). In eq. ( |6.76D D^{0) is linear in the cut-off in 
d = 3. This power divergence is absorbed in the bare parameter leading to a 



cut-off independent, universal relation eq. ( |6.78 
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The procedure is very similar in case of the spin stiffness. A possibihty to study 
the spin stiffness is to calculate the two-point function of Goldstone bosons. This 
Green-function should have a pole at A;^ = (for Hg = 0) which fixes the wave 
function renormalization factor z. For the two-point function one obtains 



-ikx 



6 



6 



-Z 



H. 



Z6Hl{x)5m{0) 
The residue of the pole gives z: 



1 Ps 



r 



po I °" ^ ^ 2d 



1 

2rf 



(6.79) 



(6.80) 



Using the Ward identity [||, [n|] 



^1/2 ' 



(6.81) 



eqs. ( |6.76| , |6l80D determine the relation between ps and 



Ps = P. 



1 - 



Ps 



(iV-2)DW(o) + 



2d 



d = 3. 



(6.82) 



Repeating the calculation for the two-point function in eq. (|6.79|) with Hg ^ 0, the 
pole position is shifted from = which determines the Goldstone boson mass. 
One obtains 



(6.83) 



In the limit Hg —>■ 0, eqs. ( |6.76| , |6!8^ , |6.83[ ) are consistent with the Ward identity |Tl 



Hs~^0 Ps 



Hg 



(6.84) 



The residue of the pole (j:'^ + (M^)^^ in the two-point function gives the wave- 
function renormalization zh, and the Ward identity [|ll| 



(M^) 



(6.85) 



gives then . The wave-function renormalization zh has the form of z in eq. ( |6.80| ) 
replacing D^^^ by D^^''"\ Eq. ([a85| ) leads to the result 



pf = Ps 



N -2 

1 + 2 M 

8npg 



(6.86) 
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Expressing eq. ( |6.83|) in terms of the renormalized parameters we get finally 

N 



in 



3 M 
ps 



(6.87) 



After restoring dimensions, eqs. ( |6.86|) and (|6.87| ) give the result quoted in eq. ( p.8|) . 

The procedure above shows that the power divergences can be consistently ab- 
sorbed in the bare parameters, the physical predictions are free of divergences as 
a — i> 0, and they satisfy the Ward identities. 



7 The effect of a weak staggered field 



Consider an euclidean box L'^^^ x Lf with "cubic" geometry, i.e. Lt/Lg =0(1). If the 
staggered field is zero, the total (staggered) magnetization moves freely in the 0{N) 
group space. Switching on an external field Hg = (iJ^, 0, . . . , 0), the distribution of 
the total staggered magnetization will be concentrated around the 0-th direction. 
For weak fields 

this distribution remains broad and there is a need to treat the corresponding mo- 
tion non-perturbatively. The external field dependence of the free energy, staggered 
magnetization and susceptibilities with O(A^) symmetry have been worked out in 



ref. W^ up to next-to-leading order. After restoring the spin wave velocity depen- 



dence we enumerated the results in eqs. ( p.21[ - pl26D . 



Increasing the box size Ls at fixed temperature, the staggered susceptibility has 
the form 



ln(L,T) + ^2.852929^ + . . .| . (7.2) 



This behaviour remains valid until Lg <^ ^{T), where C,{T) oc exp(6/T), with b = 
2-nps/{N — 2) is the nonperturbatively generated correlation length [l^, |l^. For 
infinite volume (L^ ^ HT)) the spatial size is expected to be replaced by the 
correlation length. 

These results were derived under the condition T =0(1/Ls) and, in general, the 
T 0, Lt =fixed limit does not lead to the correct zero temperature behaviour. 
The physical reason was discussed in detail in Section 6.5. As we argued there, the 
result of eq. ( p. 261) in this limit should agree with that given by the high temperature 
limit of an 0{N) rotator. We shall discuss this problem in the next section. 
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8 Low temperature limit of the finite volume stag- 
gered susceptibility 



As we discussed in Section 6, if the temperature is very low (T <^ l/-^s), the 
excitations of an 0{N) rotator dominate the partition function. Let us check first, 
how the rotator problem in eq. ( |6.47| ) is modified in the presence of a small staggered 
field Hs = {Hs, 0, . . . , 0). In the effective action a new term enters 

->[°H,-^S., (8.1) 



where is the bare staggered magnetization. Using eqs. (|6.5H6.12|) this term can 
be written as 



8.2 



Considering the first correction we can replace 7f^(x) by {l/p^){N — 1)D*{0), and 
eq. (|8.2| ) leads to a new term in the exponent of eq. ( |6.47| ): 

i J deit) exp {-^©e// (e(t) - e{t - l)f + h,e(t)| , (8.3) 



where 



,0,...,0), hs = H,Mt^fVs, Ml" = M':(^1-^^D*{0) 



The corresponding Hamilton operator reads 

1 



H 
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h.X = Ho- h,X, X 



eff 



(8.4) 



1.5) 



where we suppressed the /is-independent constant Eq{Ls) of eq. ( |6.6Up . Notice that 
this effective rotator picture is valid only in first order, i.e. including 0(1/ psLs) 
corrections. This is the order considered in this section. The staggered susceptibility 
is given by 



T 

X.(T,L.)=-^lnZ 



(8.6) 



hs=0 



where 



Z = Tr ie-^l^ 



•7) 



In taking the derivatives in eq. (|8.6|) we have to be careful since X does not commute 
with Hq. One obtains 



Xs 



Ml^f) y.— Tr / 
^ Jo 



rfAe-^^»Xoe-(i/^-^)^«Xo 
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where is the partition function with i^o? and Xq is the 0-th component of X. 
The excitations of ifo are described by the quantum number j = 0, 1, . . . 

1 



j{j + N-2) 



while the corresponding multiphcity is given by 

(j+iV-3) 



19) 



9j 



m-2)\ 



(2j + N-2). 



(8.10) 



We shall denote the set of quantum numbers which characterize the subspace with 
a given j by m. Inserting a complete set of states in eq. (^.81) we get 



where 



^0 j=o ^3 



^ |(j + l,m|Xo|j,m) 



f8.111 



(8.12) 



In deriving eq. ( ^.11 ) we used that Xq connects the state |j, m) with \ j ± 1, m) only. 
Using the relation 



1-2 - 1 



(8.13) 



J' ,m' 



and observing that the left hand side of eq. ( p.l3| ) is just aj + flj+i, we obtain 

1 (j + iV-2)! 

N j\{N~2)\ ■ ^ 
Writing eq. (|8.14|) into eq. (|8.11|) , after some algebra we get 



8.14) 



Xs 



4 2 1°*^ 



(iV-l)(iV-3) 
o^.^o 8e,ffE, + {N~l){N-3) 



(8.15) 



where t = l/{2QeffT) and U{t) satisfies the differential equation 
For = 3 eq. ( 8.15| ) simplifies since in this case U{t) = 1. 



^.16) 
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Let us consider the effective parameters JvV'J^ and 6e// of the O(A^) rotator. 
Using the definitions, eqs. ( |6.41|J6.45|J8.4| ) and replacing the bare spin stiffness and 
magnetization by their renormahzed values, eqs. ( |6.76| , |6!82D we have 



e 



PsVs 



1 + ^(Z}(0)-D*(0)) 
1 + ^^P(0)-D*(0)) + -^'' 

Ps Ps 



(8.17) 



l-d 



2d 



AoD*(0) 



where -D(O) is the lattice propagator at infinite volume and zero temperature. In 
eq. (|8.17|) we restored the lattice unit dependence to see clearly how the ultraviolet 
divergences cancel. As eqs. (A7,A14) in Appendix A show, -D(O), D*(0) and Ao-D*(0) 
are ultraviolet divergent. On the other hand, the combinations 



—a 



(d-i) 

1/2 ' 



1) + 2 



d-V 
d-2 



.l-d 



entering the effective parameters are finite. 



^.18) 



Eqs. (|8.15| - ^rT6D allows us to calculate Xs both for t <^ 1 and t ^ 1 analytically. 
(Numerically, of course, it is not difficult to perform the calculation for arbitrary 
t.) At very low temperatures, the j = term contributes only in eq. (|8.15|) and in 
Zo{t), giving 

4 



Xs{T, Ls 



N{N-1) 



TKp. <1. (8.19) 



or 



Xs{T, Ls 



MlpsV,' 



N{N-1) 
TVsPs < 1 . 
For the relevant case, c? = 3 we get 



1 + 



2N 



-a 



(d-i) 

1/2 ' 



1) + 2- 



d-l 



(8.20) 



N{N -I) 



mIpsL 



2N -3, 
A-npsLs 



-3.900265 



TLIps < 1 , = 3 . 

(8.21) 

At large temperatures TVsPs ^ 1 (but TLg <^ 1) we can use Seeley's expansion 
eq. (|6.62|) and integrate the differential equation in eq. ( ^.16|) . We obtain 



Uit) 



Zo{t) 
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■t 



-{N-3)/2 
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(A^-l)(A^-3)(5A^2_22Ar+18) 2 

t -\- . . . 



N-1 



1 - 



360 



(8.22) 



180 
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Up to leading order corrections this leads to 

XsiT, Ls) 

For d = 3 we get 

1 .LU 

1 



N 



i2ree//J 



N 



psL. 



-^ + —3.900265^+... 

47r ) 



12L,T 



(8.23) 



(8.24) 



According to eq. ( ^.21|) , in the limit T — > 0, =fixed, the staggered suscep- 
tibility becomes independent of T and is proportional to the square of the spatial 
volume. As the temperature is increased to T ^ 9 = PsK, but T <C l/i^s, the rota- 
tor becomes highly excited, while the fast Goldstone modes (whose typical energy is 
~ l/-^^s) remain frozen. In this region the susceptibility in eq. (|8.24|) should match 
the low temperature limit of eq. ( f^.26D . Indeed, using the asymptotic expressions for 
Lt ^ Ls given in Appendix B, one finds that, in the order considered in this section, 
the two expressions coincide up to 0(exp(— tt/L^T)) corrections, coming from the 
fast modes which are included in the cubic geometry. It is interesting to note that 
the susceptibility given by the effective rotator, eq. ( p.24|) reproduces the leading 
correction in eq. ( p.26| ) surprisingly well even at temperatures as high as T ~ 
(See the discussion in Appendix B.) 



9 Uniform susceptibility 



The uniform magnetic field hu in eq. ( p.l|) is coupled to the generator of the 0(3) 
rotation around the direction of h„. For general A^, will be coupled to an O(A^) 
generator Q, say Qyi generating a rotation in the 1-2 plane. We have ((5S)i = ^5*2, 
(QS)2 = —iSi, the other components of QS being zero. In the effective Lagrangean 
formulation, eq. ( |4.2| ), the uniform field enters as the time component of a con- 
stant imaginary gauge potential: 

do^ Do = do-KQ. (9.1) 

The effective Lagrangean in eq. ( |4.2| ) receives an extra contribution 

6C^ = -Kps (9oS(x)gS(x)) + ]^hlpi {QS{x)f . (9.2) 

For small /i„ it is again important to treat properly the zero modes. Let us consider 
first the "cubic" geometry. Introducing the global rotation Vt: S(a;) = r2R(a;) one 
has to average exp(— / 5Cu) over the O(A^) rotations Vt. This gives a contribution 
to the free energy 

^fu = -^As (9.3) 



hlplih I I {doR{x)R{y)){doR{y)R{x))) +0{h 



N{N-1) \ V 



y I hu=o 
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Here we have used the following relations for averaging over the group space: 

{a[b[)n = ^(ab) , (9.4) 
{{a[b', - a',b[){c[d', - c',d[))n = [(ac)(bd) - (ad)(bc)] , 

where a' = Qsl, h' = Qh, ...are rotated O(A^) vectors. The expectation value in 
eq. ( |9.3|) can be evaluated using standard steps of chiral perturbation theory. One 
obtains for the uniform susceptibility at = the result quoted in eq. (|2.29| ). 

In eq. (|2.29| ) a new shape coefficient appears defined as 

= - I d'xdlG{x)G\x) . (9.5) 

Here G{x) = G{x)+(3i{l)/ L = l/(47r|x|) + 0(x2), where G{x) is defined in eq. {^J§. 
For I = 1 ip{l) is related to /^^(l) (eq. ( |2.31| )), while for small / its asymptotic form 
is given by 

V^(/) = --i + 0(e--/n, (9.6) 



For generic values of / it should be evaluated numerically using methods discussed 
in detail in refs. ||3^, |3^. In Appendix B we give approximate analytic expressions 
for '?/'(/). 

A small staggered field Hs will fix the direction of the staggered magnetization 
provided HgAigL^ ^ T. Relative to this direction we can define the transversal and 
longitudinal uniform susceptibilities, XuC^^^s) and xi{T,Ls), respectively.|^ The 
susceptibility Xu(T,Ls) given in eq. ( |2.29| ) corresponds to Hg = 0, and is given by 
the "angular" average over the O(A^) group space 



(iV-l)x-(T,L.). (^-r-^' x!l(T,L. 



(9.7) 



In the infinite volume limit Xu{T) and Xu{T) are infrared divergent in chiral pertur- 
bation theory. On the other hand, the O(A^) invariant average Xu{T) remains finite 
and it is given by eq. (|2.32| ). 

Small temperatures T ^ l/Lg should be considered again separately. The 
Hamiltonian for the effective O(A^) rotator in this case is 

H=-^V-KU, (9.8) 



^^For the general 0{N) case "parallel" is defined by the 0{N) generators under which the 
staggered magnetization is invariant. 
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where we denoted for simplicity by the O(A^) generator coupled to With the 
notations of Section 8 we obtain 

For TVsPs ^ 1 (but T ^ ^/Lg) Seeley's expansion gives 

^ _ 1 M ^ 

, (9.10) 



Xu{T,Ls) = —ps < 1 + 



in agreement with the low temperature behaviour obtained from eq. ( |2.29D . 

For very low temperatures, TVgPs <^ 1, only the lowest lying rotator states 
contribute, leading to an exponentially small susceptibility given by eq. ( |2.33| ). 
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A Momentum sums and integrals on the lattice 

It is useful to introduce the notations 



= e'^^ - 1 , k*^ = e-""^ - 1 , d{k) =4J2 sin' ^ 



- (27r: 



id 



r=E> e;= e , (A.i) 



k k k 

fcO/O fc7^(0,...,0) 



1 



d{k) 

With the help of the identities 



one obtains 



(kfj, + k*A = kjx,*^ = 4 sin^ , (no sum over fx) , (A. 3) 



= -^E*^^^^(^ = -2A^^*(0). (no sum over /i) (A.5) 
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In the following we shall discuss the volume dependence of D*{0) and A^D*(0) in a 
euclidean box L'^~^ x Lt in the ^ oo limit. 

The sum in D*{0) can be written as 



lim 



E 



d{k) [ V Y + rf(A;0, A; = 0) + 



(A.6) 

We divide the sums in eq. ( [A .61 ) into the infinite volume contribution and the rest, 
which is volume dependent and remains finite as the lattice unit a —>■ 0: 



1 W_l 



lim 



1 



ko 4 sin^ ^ + M2 



(d=i) 



(A.7) 



Since we are not interested in cut-off effects which go to zero as the cut-off goes 
to infinity, we can evaluate gi^^ in the a — *• continuum limit. This problem is 
discussed in detail in the Appendix B of ref. [ITll and we quote the result only: 



a? 



1 



47rLf-2 
0. 



a 



(d-i) 

1/2 



(1) + 



Ait 



2ML, 



- 2- 



d 



d-2 



(A. 



In the M — > limit the first integral in eq. ( [A.7[ ) is finite while the second one gives 
1/2M. We obtain 



+ 



kd{k) A-kL^ 



a 



(d-i) 

1/2 > 



l)-2 



d 



d-2 



(A.9) 



The first term is a volume independent but cut-off dependent constant, while the 
number afl^^\l) can be found in Table 1 in ref. [pH] . 



The steps are similar for A^D*(0): 

A.^*(o) = ^E* 



d{k) 



(A. 10) 



lim 



k,. 



m"Zo I V Y d{k) + M2 



Vt71~E 



K ^ d{k^, k = 0) + M2 



The second term in eq. ( |A.10| ) gives +6^o/2Vs. In the first term we separate the 
infinite volume contribution: 



1 



k,^ 



1 



M™o V^d{k) + M^ 2Jk d{k) 



(A.ll) 



where 



lim - 

Af2^0 2 



1 + k*^ 



V d{k) + M2 Jk d{k) + M2 



(A.12) 
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It is easy to show that 



d- 1 



Qo, z = l,...,d-l 



(A.13) 



Indeed, for — »• oo we can write 



1 

2 I K 



2 K 



- rfA;o-4sin2|i 
dko Y.i 4 sin 



(A. 14) 



2 ^ 
2 



27r d{k) 



Since we have a box ^ x Lt, Qi is independent of i and eq. (|A.13|) follows. The 
first term in eq. ( [A. Ill ) gives 



1 1 



2 Jk d{k) 



2d a'^-i 



(A.15) 



where on the right hand side we restored the cut-off dependence explicitly. We have, 
therefore 



AoI^*(0) 
A,D*(0) 



11 1 „ 
■2d^ + 2V: + ^° 



(A.16) 



1 1 



2c? a 



d-l 



d-1 



Qo. 



In the expression for the ground state energy eq. ( 6.57 ) the combination 



1 1 



2d a 



d-l 



2K 



(A.17) 



1 1 



2d a-^-i d - 1 
enters. Eq. ( |A.17|) leads to the result 



1 



Qo -77 



1 1 



K V 2d a" 



A,D*iO)A,D*iO) - -AoZ^*(0) = -- 



1 . 11 d 1 



2d-2 



+ 



d- 1 



(A.18) 



As expected, there is no mixing between divergent cut-off powers and volume de- 
pendence in eq. ( [A.18|) . Since Qo is finite in the a — limit, we can write 



lim 



°° dk^ kl 



27r B + fc2 



(A. 19) 



1/2 



1 



40F^ 



id-l) 
'-1/2 
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where, in the last step, we used the notation of ref. [0. Since 
restoring the dimensions we get 



(d-i) 

/2 



after 



Qo = 



a 



and can be neglected in eq. ( A.18| ) in our order. We can write therefore 



A^D*(0)A^Z}*(0) - -AoD*(0) 



1 1 



(A.20) 



(A.21) 



The divergent first term contributes to the non-universal constant of the ground 
state energy density only. 



B Asymptotic expressions 



For the "slab geometry", / <C 1 (or Lt ^ Lg) and for the "cylinder geometry", 
/ ^ 1 {Lt ^ Lg) one can derive simple asymptotic expressions for the functions 
ar{l) = af\l)^ which are related by eq. (|2.11|) to the functions n = 0,1,2 

appearing in the finite size corrections. One can show that for / < 1 



af\l) 



I 



4r-6 



„(1) 



1 



+ O (e--/'') , / < 1 . 



+ 1 



-2r 



r(2r - 3) 

For / > 1 one has the approximation: 

1 



I 



4r 



a 



r(2r - 1) 



3-2r 



a 



(2) 

-r+3/2 



(1) 



r(2r - 3) 



blT + 0(e-"). l>l 



r{r — 1) 



(B.l) 



(2r - l)(2r - 3) 



(B.2) 



Using these expressions and the equality a^\l) = ttrf/2-p(-'-) relate the 

shape coefficients to those corresponding to a symmetric box in d dimensions, listed 
in Table 1 of ref . 



In particular we have for / < 1 



■ 0.382626 + 2 In / + 1.054689 + Sp^il) 
^ [2.852929 + 61n/] + y(5/?f (/) , 



{Any 



0.610644 



/6 



+ 1.047198 



+ J^S(3^{1) 



(B.3) 



40 



n 




K 




K 





8.85264 


-0.959427 


8.62200 


-0.608985 


1 


-0.66610 


-0.321512 


-0.671988 


0.553185 


2 


0.058438 


1.35468 


0.061051 


1.25040 



Table 1: Values of the parameters in eq. ( B.5 ). 

For / > 1 we get 

/3o(/) = 1.437745 -4 In/ + 5/^0+ (0, 
1,3 3.900265 1 ,^+,,, 

The correction terms S(3^{1) are exponentially small for / ^ 1 and / ^ 1, respec- 
tively. They can be parametrized in a simple form 

a-^^e-2-/'\ /<1, 
6(3:{l)=a+l''-e-'-^\ 1>1. (B.5) 

The exponential factors are expected on physical grounds, while the parameters a J, 
can be determined from the known values of /3n(l) and the fact that /3^(1) = 0. 
They are summarized in Table 1. It is remarkable that this simple parametrization 
gives a relative error smaller than 10~^ in the whole range of / values. 

For the shape coefficients ■?/'(/) defined in eq. ( |9.5D one can also give approximate 
analytic expressions 

^(0 = -^ + 2.57e-2-/'% /<1, (B.6) 

^(/) = + 3^900265^5 _ ^^^^ ^2 _ ^ 
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